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MOTION  OF  A  ROCKET  DURING  THE  BURNING  PERIOD 


Abstract 


— In  thif report  the  general  equations  of  motion  of  a  rocket  are 
developed.  These  equations  are  then  applied  to  a  rocket  which  is 
subjected  to  an  eccentric  propelling  force,  one  of  the  principal 
causes  of  inaccuracy  in  rocket  firings.  The  equations  are  first 
specialized  for  the  vacuum  case  of  the  non-rotating  rocket.  The 
results  for  this  case  are  compared  with  those  for  a  rotating  rocket 
moving  in  a  vacuum.  A  formula  for  the  angular  velocity  necessary  to 
reduce  dispersion  is  given.  Finally,  aerodynamic  forces  are  con  - 
sidered  and  expressions  for  the  yaw,  angle  of  deviation  and  deflection 
are  set  •’own.A 


1.  Introduc  tion 

The  equations  of  motion  which  govern  the  flight 
of  a  symmetric  rocket  will  be  developed  in  a  general  way 
in  order  that  they  may  be  used  for  reference  purposes,  and 
used  in  the  analysis  of  problems  associated  both  with  the 
burning  period  and  the  post-burning  period. 

After  the  equations  have  been  set  up,  this  report  will  be 
primarily  concerned  with  the  effect  of  an  imperfect  alignment 
of  the  propelling  jet.  The  eccentricity  of  the  propelling 
force  is  cne  of  the  main  causes  of  inaccuracy  in  the  firing 
of  rocket?.. 


2. 


Notation 


In  setting  up  the  equations  of  motion  for  a  rocket 
the  following  symbols  and  conventions  will  be  used: 

—l*  Co>  Co  =  unit  vectors  defining  a  right-handed  Newtonian 
reference  ^frame.  The  vector  £,  points  in  the  direction  in 
which  the  rocket  is  launched.  x 

i,  At  k  =  unit  vectors  defining  a  right-handed  moving 
reference  system  whose  origin  0  is  at  the  center  of  mass  of 
the  rocket.  The  vector  .i  coincides  with  the  longitudinal  axis 
of  the  rocket.  The  vectors  j_  and  k  are  temporarily  left 
undefined.  Unless  there  is  a  statement  to  the  contrary,  all 
vectors  are  consider? d  as  resolved  with  respect  to  the  moving 
axes . 

v  =  (v^,  v2,  v,,)  =  velocity  of  0. 

o)  =  (co^,  (0^,  co3)  =  angular  velocity  of  rocket. 

J!_  -  Q2>  £3)  =  angular  velocity  of  moving  axes. 

3?  -  (F1,  F2,  Fj)  =  resultant  propelling  force  at  0. 

P  =  (P^,  P2,  P^)  =  resultant  propelling  moment  about  0. 

f  =  (f^,  f2,  f^)  =  resultant  aerodynamic  force  at  0. 

P  =  (P^>  P2,  P3)  =  resultant  aerodynamic  moment  about  0. 

_g  =  (g^g^*  s3^  =  acceleration  due  to  gravity. 


m  =  mass  of  rocket. 

A  =  axial  moment  of  inertia  of  rocket. 

B  =  transverse  moment  of  inertia  of  rocket. 


N  =  axial  spin  of  rocket  relative  to  moving  axes. 

_h  =  Ja(N+  B  Q2,  BC3]  =  angular  momentum  of  rocket 
about  0.  Since  R2=co2,  =  <o3 ,  the  angular  momentum  can 

be  written  *' 


h  =  (Aa>^,  Ba>2,  Bco^) . 
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3. 

are 


General  Equations  of  Motion 

In  vector  form,  the  fundamental  equations  of  motion 


dt 


Written  out  at  length  in  scalar  form  these  equations  yield 
(m'v1)  +  ra(a)2v3  -  v2co3)  =  F^  +  f^  +  mg1 

(m'v2)  +  m(a>3v1  «  v^)  =  F2  +  f2  +  mg2 

t 

(m'v3)  +  raCQ^  -  v^)  =  F3  +  f3  +  mg3 

(A  i1)  -  P2  +  px 

(Bu>2)  +  Aaj^  -  Bco^  “  P2  +  P2 

(Ba>3)  +  BQ^  -  Aax^cOg  =  +  p3 

where  the  dot  denotes  differentiation  with  respedt  to  the 
time  t.  The  above  equations  can  be  expressed  more  conven¬ 
iently  and  more  compactly  by  introducing  the  complex  quantities 

vc  =  v2  +  iv3  °>c  *  “>2  +  ic°3 

Fc  =  F2  +  iF3  f c  =  f2  +  if 3  Sc  =  .g2  +  ig3 

Pc  =  P2  +  ip3  Pc  ”  P2  +  ip3* 

In  terms  of  these  quantities  the  equations  of  motion  reduce  to 
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(m'v^  +  m(a)2v3  “  VV  =  Fi  +  fi  + 


(ra#vc)  +  "  vi«c)  =  Fc  +  fc  +  ragc 


(AcDj^)  =  F2  + 


(^c)  -  i  (A^  -  BQx)  a>c  =  Pc  +  Pc- 


For  certain  parts  of  the  development  which  follows  it  is 
necessary  to  know  the  moving  axes  components  of  a  vector  r 
which  is  fixed  with  respect  to  the  inertial  frame.  If  this 
vector  is  written 

-r  =  ri  i  +  V  +  r3^ 

the  components  can  be  found  by  noting  that  r  =  0  which  leads  to 
*1  +  "2r3  '  “3r2  =  0 

r2  +  “2rl  -  ^1*2  =  0 


r3  +  Qir2  ~  Vl  =  0 


or,  if  rc  =  r2  +  ir^  is  introduced,  to 

*1  +  “2r3  -  V2  =  0 
■  *c  +  1(Qlro  -  rl“c)  *  °- 


These  equations  associated  with  the  proper  initial  conditions 
for  the  particular  vector  r  determine  its  moving  axes  com¬ 
ponents  r^,  r2>  r^. 

It  will  be  interesting  later  to  find  the  yaw  5,  (the  angle 
between  i  and  v) ;  the  angle  of  deviation  9,  (the  angle  between 
v  and  cji  and  the  angle  f(p  between  _i  and  c^.  These  angles 

can  be  found  provided  v  and  Q  are  known.  The  yaw  is  given  by 
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tan  6  = 


i  x  v 

V2  k  -  V3I  ' 

\I4  *  v3 

vc 

i  •  v 

V1 

V1 

V1 

If  r  = 


the  angle  of  deviation  is  given  by 

tan  9  = 


P-1  x  Z 


% '  v 


Wlr2  ~  r2vli2trlY3-r3vl)2'f(r2v3~I'3v2) 


rlvl  +  r2v2  +  r3v3 


which  for  sufficiently  small  r5,  r.,,  v9,  v^  can  be  replaced  by 

* - * -  ^  J  ..u'-iU-J 


tan  0' 


vlvl 

1  Vl 

v2  *2 

vc  rc 

V1  rl  + 

V1  rl 

V1  rl 

(rlv2-r2vl> ,  (rlv3-r3 


rlvl 


1 


The  angle  (p  is  given  by 
tan  <p  = 


i  X  £1 

]/4  +  r^ 

rc 

i  .  c^ 

rl 

rl 

For  the  usual  burning  period  the  above  angles  are  generally 
small;  the  tangents  of  these  angles  can  then  of  course  be 
replaced  by  the  angles  themselves. 

4 .  Aerodynamic  Forces 

The  analysis  of  the  aerodynamic  forces  acting  on 
the  racket  will  be  based  on  the  assumptions  and  methods  in¬ 
troduced  by  Nielsen  and  Synge  in  their  very  important  paper.* 

In  accordance  with  these  authors,  it  is  assumed  that  each 
component  of  f  is  a  function  of  p,  c,  v1,  v2,  v^,  co-^,  o>2,  co^; 

where  p  is  the  density  of  the  atmosphere  and  c  is  the  velocity 
of  sound.  If  v2,  Vy  o>2,  o)j  are  small  the  components  f2  and 

f^  can  be  replaced  by  the  approximations 


*  K.  L.  Nielsen  and  J.  L.  Synge,  On  the  Motion  of  a  Spinning  Shell, 
reproduced  by  the  Ballistic  Research  Laboratory  by  permission  of  the 
National  Research  Council  of  Canada,  1943* 


***** 
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f2  =  alV2  +  blV3  +  cl®2  +  dla>3 


=  a2v2  +  t»2v3  +  c2a>2  f  d2a>3 

in  which  constant  terms  do  not  appear  since  if  v2  =  v^  =  o>2  =  ©3  ~  0 
then  f2  and  f^  must  vanish.  The  coefficients  in  the  approxi¬ 
mations  are  functions  of  p,  c,  a^.  The  cross  force 

f  =  f0  +  if,  is  thus 
c  2  3 

fc  =  (al  +  la2^  v2  +  (bl  +  ib2^  v3  +  (c1+ic2)(D2+(d1+id2)co3 
or  using  vc  =  v2  +  iv3>  “c  =  ^2  +  ^3 


where  the  coefficients  are  complex  functions  of  p,  c,  v1#  co-^. 

'  •  It  is  assumed  now  that  the  rocket  is  symmetric  in  the 
sense  that  a  rotation  of  the  rocket  through  the  angle 
^  =  2ir/n,  where  u  is  an  integer  greater  than  2,  about  its 
longitudinal  axis,  restores  the  rocket  to  its  original  position. 
For  such  a  rocket  the  coefficients  of  the  conjugate  terms  in  f 
are  zero.  Synge  and  Biot  have  shown  this  in  the  following 
way.  For  any  given  motion  we  can  consider  an  alternate  motion 
defined  by  velocity  vectors  which  result  from  rotating  rigidly 
the  velocity  vectors  v  and  to  through  the  angle,  ip  ,  about  the 
axis  of  the  rocket.  If  this  alternate  motion  is  considered  as. 
referred  to  the .same  axis  used,  to  describe  the  actual  motion, 
then  v(;e1y,  a^e1*  are  the  new  velocities.  From  the  symmetry 

of  the  rocket  the  new  cross  force  must  be  f  e*^  and  according 

iv  c  ■ 

to  the  assumptions  f  e  ’  must  be  equal  to 


=  °^vceil*1  +  ^  + ^  *  W1’  • 


If  the  original  expression  for  f  is  substituted  in  this 
equation,  it  is  found  that 


2i  0 

.h  0 


2ibsin  y  =  0 


h. l  =  0 


and  hence 


fc  =  alvc  +  ^1®C 


By  the  same  analysis 


»c  =  «2vc  +  W 

The  cross  force  and  the  cross  moment  can  be  written  in 
the  more  explicit  forms 

fc  =  [^l(vi>0)i)+  ie2(vi»a)l)Jvc  +  [93^vl,a)l)+i94(vi^l) 
Pc  =  J95(vi,a)i)+i06(vi>®i)]vc  +  e7(vi*a>i)+108(vl,a)l)| 


03 


CO 


CO 


It  .  is  also  assumed  that  f-^  and  are  functions  of  p, 
c,  vlf  v2,  V3J  co^,  co2,  co3;  but  that  the  effects  of  v2,  v3; 

2,  co3  on  f^  and  p-^  can  be  entirely  neglected.  Consider  then 

f^  and  p-^  as  functions  of  p  ,  c,  v^,  co^only.  From  the  dimensional 
standpoint,  it  is  convenient  to  define  quantities  K^,  and  KA 
by  means  of  the  equations 

fl  ■  PT1Z  ^DA 


Pi  =  pv 


where  d  is  the  diameter  of  the  rocket  and  the  K’s  are  dimensionless 
functions,  f,  is  the  axial  drag  and  p^  is  the  spin  decelerating 
moment .  x  x  . 

The  9*s  v/hich  appear  in  f„  and  p„  can  be  interpreted  as 
follows.  Suppose 


co^  =  co2  =  co3  =  0,  then 
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i 

f 


(v2+iv3) 

(v2  +  iv3). 

If  in  addition  =  0,  then  f 2  =  0i(vi>°)v2i  *3  =  e(vi>0fv2> 

P2  -  03(v1,O)y2,  ;  p3  =  (^(vp  0)?£  j  but  for  this  case  the 

motion  is  planar  and  f3  and  p2  must  be  zero.  To  meet  these 
requirements  set  Op  =  co^  “1)  J  ®5  =  a>i®5*(vi>  cd^)  . 

Then  for  a)^  =  a>2  =  a>3  =  0,  v^  *  0, 


f0 

=  ®l(vl* 

vc 

?c 

=  e6(Vl,  0) 

Vc 

Dimensional  considerations  then  lead  to 

A  =  p  Tid2Ku  e6  =  p  vi  43  km- 

The  quantity  01vc  is  the  cross  force  due  to  cross  velocity, 

on  the  normal  force,  and  it  has  the  same  orientation  as  the 
cross  velocity.  The  quantity  &£VC  is  the  cross  torque  due 

to  cross  velocity,  or  the  restoring  torque. 

Suppose  now  that  v1  =  Q;a>2  =  a>3  =  0.  Then 

fc  =  *l>|vc| 

Pc  =  Wjl ©|  (0,  <o1)|  vc| 

and  from  dimensional  considerations 

02  =.pa>id3KF  0$ 


The  quantity  02v  i5_the-Magnugcross  force  due  to  cross 

"velocity,  and  the  “quantity  G^vi*  is  the  Magnus  cross  torque  due 
to  cross  velocity.,  .  I'  ■?’  c 


/ 

If  v2  =  v^  =  Oj  co1  »  0,  then 


Pc  =[©7(v10)  +  ie8(vi,0)j  (o)2  +  ico3) 

and  if  in  addition  a>2  =  0,  then  f2  =“9^(v1>  0)u>3;  f^  *  63(v1,0)a>3J 

P2  =-9g(v1>  0)a>3;  P3  =  0y(vi,O)o3J  but  for  this  case  the  motion 

is  planar.  Therefore  f^  and  p2  must  be  zero. To  meet  these 
requirements  set  0^  =  co^G^*(v^,  a^) $  0g  =  co^0g*(v^,co^) .  Then 

for  v2  =  v^  =  0,  co-j^  -  0, 

ro|  =  tyV^c 

l 

1 5>C  |  =  ®7(tx>0>  “c 

Dimensional  theory  then  leads  to 

e4  =  Pvid3ks  e7  =  PTia\- 

The  quantity  coc  is  the  cross  force  due  to  cross  spin  and 

Ine  quantity  0„co  is  the  cross  torque  due  to  cross  spin. 

■  c 

Suppose  now  that  v2  =  v^  =  0,  v^  =  0.  Then 

|fc|  =  ule3*(0'  “1)  |“C 

|  Pc  |  =  “i08*(°'  “l>  |“c 

and  again  from  dimensional  considerations 

93  =  6g  =  pOjd5!^ 


The  quantity  03cuc  is  the  Magnus  cross  force  due  to  cross  spin, 
and  the  quantity  0goac  is  the  Magnus  cross  torque  due  to  cross  spin. 


5 .  _  Vacuum.  ^&se  — , ; ' 

It  is  instructive  to  study  first  the  trajectory  of 
a  non-rotating  roc.cet  which  is  subjected  tp  a  propelling 
force  F(t)  which  acts  parallel  to  the  longitudinal  axis  of 
the  rocket  at  a  distance  of  e  units  from  that  axis.  In 
this  section  the  effect  of  such  an  imperfect  alignment  of  the 
propelling  jet  on  a  rocket  moving  in  a  vacuum  will  be  in¬ 
vestigated. 

Let  the  moving  axes  be  fixed  in  the  rocket  and  take  the 
Jc  axis  so  that  it  coincides  in  direction  with  the  propelling 
movement  vector.  Then  f  =  O  j  £  =  0;  F^  =  F^(t)j  F2  =  Oj 

F3  =  0;  =  0;  ?2  =  °i  and  P3  =  For  the  particular  case 

under  consideration,  it  is  convenient  to  introduce  £  =  vi+iv2> 

=  gjL  +  ig2  and  r)  =  co2  +  ico^.  Taking  into  account  the  fact 

that  the  mass  of  the  rocket  varies  with  the  time,  the  equations 
of  motion  are 

(Aa>,)  =0 

(Bn)  -  i(A-B)oo1n  =  IF^ 

(m*£  )+  imco^  £  +  mv3(<D2-i<D^)  =  Fl  +  mcf4 
(mv3)+  m(co1v2  -  v^)  =  mg^ 

If  time  is  measured  from  the  instant  at  which  the  rocket 

emerges  from  the  launching  tube,  the  initial  velocities  are 

v  =  v  i:  co  =  0  for  t  =  0. 

—  Or—  — 

The  equations  for  the  components  of  £  are 
cf  +  ico^  <7  +  g^  (co2  -  ico^)  -  0 

g3  +  -  a)2Sl  =  0 


The  initial  values  of  the  components  of  £  depend  upon  the 
choice  of  the  inertial  frame.  This  frame  is  chosen  so  that 
c^  coincides  with  the  axis  of  the  launching  tube  and  c.^  is 

horizontal.  Let  (x  be  the  angle  of  elevation  of  the  tube. 
The  components  of  £  referred  to  the  in  .  .  frame  are  then 
(-gsinji  ,  -  g  cosp  ,  0).  Suppose  that  the  >lane  determined 
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by  i  and  F(t)  at  t  =  0  makes  an  angle  v  with  the  plane 
determined  by  and  With  respect  to  the  moving  ax&s, 

the  initial  components  of  £  are  then  (-g  sin  p,  -  g  cos  p  cos  v, 
g  cos  p  sinv  ) . 


From  the  first  equation  of  motion  cu^  =  0  and  hence  from 
the  second  ico^)  =  i  J  F^  dt  =  i  h(t)  so  that  co^  =  0; 

'o 

=  h/B.  The  components  of  g  are  given  by 
g^  =  cons 1 1. =  -  g  cos  p sin  v 

and 


cr  +  i  oj^ct  “0 


from  which 


'  -/: 
=  gl  +  i?2  =  Cle 


co_dt  •  / .  \ 

3  = 


where  s(t)  =  J  co^dt  =  and  C-^  =  -gsin  p-  igcosp  cos  v. 


The  remaining  equations  reruce  to 

-  -R  cos  u  sin  v 


V3  = 


m 


/; 


dt 


and 


(m£  )  +  i  m  cu^  5  =Fi+  nC^e"^s^^ 


from  which 


m  £  =  m(v1  +  iv2)  =  j  F1e 

/o 


[s(T)-s(t)J 


a 


C-,  m  dt  +  m  v^ 
1  o  o 

0 


dT 


e 


-i  s (t) 


where  mQ  is  the  initial  mass  and  vQ  is  the  initial  velocity. 

If  r  is  a  vector  which  is  fixed  with  respect  to  the 
inertial  frame,  the  components  r^,  r^>  r^  of  this  vector  with 

respect  to  the  moving  axes  are  given  by  the  equations 
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®3r2  =  0 


« 

rl  " 

*2  +  Vl  =  0 


If  C 


or 


from  which  v-^,  v 2  and  are  easily  determined.  Applying 
these  equations  to  £^,£2  anc*  —3  we 


Initial  components 


<y  * 

-1  • 

«'»  • 

“2  * 

(1> 

(0, 

0,  0) 

COS  V  , 

sin  v  ) 

K  =  e15^ 

Z,  -  -i  cos  v  eis 

H  II  ■ 

u  u 

0 

sin  v 

£3  : 

(0, 

-sinv, 

cos  v  ) . 

Z  =  i  sinv  eis 

r3  = 

cos  V 

If  £  is  a  unit  vector,  the  component  of  v  along  r  is 
v  .  r  =  v-jT^  +  v2r2  +  +  v^.  Therefore  the 

components  of  v  with  respect  to  the  inertial  frame  are: 


=  r^  -  ir2  is  introduced,  these  equations  become:; 


Mfl)5(  =0 


r3  ■« 


C  =  C2eis(t) 


r^  =  cons»t. 


V1  =  i cos[s(Tj]  dT  -  JLMgULj  ^ 


v2 =  ttv/t  sin[s(T!3 dT  -  £-r£tlf“ 

=  -  Sify- Jh  slnf(T^  dT' 


F1  sin  s(T)  dT 


Writing  M(t)  =  /  —  the  integrals  of  the  above  equations  can 
be  written 

t  t 

X  =  J  M(t)  -  M(T)j  Fx  cos[s(T)JdT  -  g  sin  p  /  M(t)-M(T)j  ra(T)dT 

^  /> 


+  m0v0  M(t) 


Y  =  cosv 


;^(t)-M(T)j  F1  sin{s(T)J  dT  -  g  cos  p^M(t)-M(T)|  m(T)dT 


Z  =  -  sin  v 


Thus  the  restricted  problem  of  this  section  is  solved  in  terms 
of  quadratures. 

The  deflection,  and  the  angle  which  the  axis  of  the  rocket 
makes  with  the  line  of  departure  (the  £,  axis)  are  given 
respectively  by  the  equations 


The  angle  of  deviation,  and  the  yaw  are  given  respectively  by 
the  equations 

tan  9  =  \/v|  +  V^/V-l  ;  tan  b  =  ^Vj  +  V|/Vr 

All  of  the  above  quantities  except  <p  involve  gravity  terms. 
These  terms,  however,  can  generally  be  neglected. 

If  ra,  F,  B,  and  e  are  regarded  as  constants,  and  if 
gravity  is  neglected  the  deflection  d  is 
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the  angle  <p  Is 


and  the  angle  of  deviation  is  r 


tan  y  = 


=  11 
m  l 


C  F-jEt- 

sin  —4^ —  dt. 


2B 


Approximate  expressions  for  d  and  9  are 

F18et^' 

d  =  24  m  B 

Fj&t3 

0  =  6m  BV 


The  yaw  is  the  difference  <p-9.  Thus  it  is  seen  that  the 
quantities  d,  q>,  0  fe  are  approximately  directly  proportional 

to  the  eccentricity  factor  e. 

6.  Rotating  Rocket  Vacuum  Case 

If  a  rocket  is  caused  to  rotate  about  its  longitudinal 
axis  with  an  angular  velocity  co^,.it  is  not  unreasonable  to 

suppose  that  will  tend  to  overcome  any  eccentric  action 

of  the  propelling  forces.  In  order  to  investigate  this  suspected 
stabilizing  effect  of  co^,  consider  now  the  motion  of  a  rotating 

rotket.  For  the  time  being  suppose  that  the  motion  takes 
place  in  a  vacuum,  and  suppose  that  the  effect  of  gravity 
can  be  neglected.  Let  the  propelling  force  F  and  the  pro¬ 
pelling  moment  P  be  arbitrary  functions  of  t.  Let  the  moving 
axes  be  fixed  in  the  rocket,  and  take  k  on  the  plane  deter¬ 
mined  by  P  and  j..  Choose  the  inertial  frame  so  that  c  coincides 
with  the  axis  of  the  launching  tube,  and  c*  is  horizontal.- 
Let  the  rocket  emerge  from  the  tube  where  t  =  0.  and  assume 

that  for  t  =  0,  v  =  v  i,  co  =  (pi. 

o~  —  o— 

The  equations  of  motion  are 

(A^)  =  Px 
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(B^c)  -  iCA-Bjo)^  =  iP3 


(mv^)  +  m  (o>2v3  “  V2C03)  s  Fi 


(mVc)  +  icu^  mvc  -  imv^ 


=  i’  • 


From  the  first  equation 


“i  ■  I  f  pi 


dt  +  CD  . 
o 


Regarding  Bar  as  dependent  variable  in  the  second  equation, 

V 

the  solution  of  this  equation  is  easily  found  to  be 

,dt 
4  J 
e 


"c  =  B  e 


-i  j[  kco^dt  kco^c 


P3  dt 


where  k  =  1  -  A/B.  Since  F., 


and  are  small  in  comparison 


with  F1  and  respectively,  the  quantity  c^v^  -  v^co^  can  be 
neglected  in  the  third  equation  which  then  yields 


'i  -  £  fh 


dt  +  v  . 
o 


The  solution  of  the  fourth  equation  can  now  be  found.  It  is 
-ijoj-^dt  r*  i/co^dt  ijeo^dt  ^  ijeo^d l 

A 


v  =  —  e 
c  m 


F- dt  +  ~  e  * 
c  m 


mv,a)  dt 
1  c 


The  moving  axes  components  of  a  vector  r  which  is  fixed 
with  respect  to  the  inertial  frame  are  given  by 


rl  +  c02r3  "  t03r2  =  0 


rc  +  1,J,lrc  -  lrl“c  = 


If  r  =  and  a^r^  -  (03r2  ne§lected, 
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%  «**«- 


I 


=  cons' t» 


■t  t 


7s  =  i  ® 
rl 


-ij"  c-^dt  r  co^dt 


»cdt. 


Expressions  Tor  the  yaw,  the  angle  of  deviation  and  the 
deflection  can  now  be  set  down.  The  yaw  b  is  determined  by 


tanb  = 


1 

mvn 


ft,  l/to] 

6  Fcdt  -J  —  e  *  J  P3e  '  ■ 


dt 


dt  dt 


and  the  angle  of  deviation  is  determined  by 


tan  0  = 


Ic  lc 
V1  rl 


1 

DV-, 


t 

ij^dl 

e  F  dT 

v  * 


rt 

+  r 

f0 

f 


j  ~f  S 

-  if(l-k)oo,dT  T  i/kco,as 

[»(t)v1(t)-«(T)v1(T)7_  4  l_'  „ 

-  e  y  e  p: 


dsdT 


The  component  of  v  perpendicular  to  c.-,  is  v,  tan  0.  The 

/'t  x  ± 

deflection  d  is  then  J  v^  tan  0  dt. 

Suppose  that  m,  B,  Fj,  e  are  constants;  ana  suppose  Fc=0, 
©1=0,  P^=  Fj^e.  Then 


V2  p  / 

*1 


(t-T) 


ds  dT  = 


Fjet3 
6m  13 


and 


F?  et^ 
d  =  — - 


24mB 


-16- 


•  •# 


These  results  check  with  those  found  in  the  previous  section. 


P3 


Suppose  now  that  m,  B,  A,  F-,  >  e  are  constants;  F  =0, 

,  *  ^  c 

F^e,  but  that  0.  Then  if  J  co^dt  =  s(t) 


p2  g 

vi  tan  9  =  5T 


■t  /T 

i(l-k)s(T)  (  iks(x) 

(t-T)e  |  e  dx  dT 


Che  integral  which  appears  above  can  be  appraised  •  as  follows, 
integration  of  the-  inner  integral  by  parts  leads  to 


v^tanQ 


Ffe 

51“ 


(t-T)e 

ikcon 


is(T)  n  . 

dT  -  tl-Al 


i(l-k)c(T) 


»  lKO>0 


dT 


4iir  f 

*  rs 


i(l-k)s(T) 


iks(x) 

/  sf  X  dx  dT 

p  1 


Noting  that  the  absolute  magnitude  of  a  sum  is  equal  to  or  less 
than  the  sum  of  the  absolute  magnitudes  of  the  terms  of  the  sum, 


<  Ffe 

V1  tan  9  =  5F" 


-l 


(t-T)dT 

kto^ 


♦4 


(t-T)dT 

kc> 

o 


(t-T) 


7 


co; 


dx  dT 


Integration  then  yields 


<  Ffe  t2 

v,  tan  0  =  — =n - 

1  mBko> 
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—~~~r  i  ruin . . 


If  this  result  .is  compared  with  the  corresponding 
result  for  the  case  of  no  rotation,  viz.,  v^tan  9  =  F£et3/6mB, 

it  can  be  easily  seen  that  if  the  angular  velocity  is  to  be 
effective  in  reducing  the  angle  of  deviation  then 

1  <  t 

k©  ”  T> 
o 


or 


It  should  be  noted  that  if  the  above  procedure  for  the 
appraisal  of  v1  tan  9  is  applied  to  the  expression  (a)  a 

better  estimate  can  be  obtained,  namely  ; 


< 

v^  tan  0  = 


t2  .  2t  .  Xt2 

2k©^  k©2  k^co3  * 

0  0  o 


7.  Rotating  Rocket  Subjected  to  Aerodynamic  Forces. 

The  motion  of  a  rocket  during  the  burning  period 
will  now  be  investigated  under  less  restrictive  conditions 
than  have  hitherto  been  imposed.  Three  of  the  aerodynamic 
forces  will  be  admitted  into  the  equations  of  motion  of  a 
spinning  rocket,  namely,  the  drag  =  pd2KDv2  =  J^v2, 

the  lift  =  pd2KjV2  sinb  =  J^v2  sin5,  and  the  righting  moment  = 

pd^jV2  sin  ~  JMv2sin5,  where  5  is  the  yaw.  Only  these  forces 

are  admitted  because  they  are  the  most  important.  It  is  to 
be  noted,  hov/ever,  that  even  if  all  the  aerodynamic  forces 
are  admitted,  the  equations  of  motion  can  be  solved  in  the 
manner  discussed  below. 

Let  the  inertial  frame'  be  chosen  as  above.  Let  the 
moving  axes  be  fixed  in  the  rocket  with  _i  coinciding  with 
the  axis  of  the  rocket  and  k  in  the  plane  determined  by  i. 
and  P.  Suppose  that  the  rocket  emerges  from  the  tube  when 
t  =  0;  and  that  for  t  =  0  v  =  v  i,  to  =  ©  i. 

—  o-  Or- 

Using  the  terminology  introduced  earlier,  the  axial 
drag  f-^  is  equal  to 

**1  =  ^Lv2  s*n2k  “  Jj}vSccs  b 
=  J^v2  tan2£>  -  JDv£  sin 6 
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gl  +  “2g3  -  “3g2  =  0 
4  +  l“lgc  '  lgl“c  =  °' 


No  attempt  will  be  made  here  to  solve  the  above  equations 
of  motion  in  a  general  manner.  Instead,  the  analysis  will  be 
confined  to  the  simplified  equations  which  result  from 
assuming  a^v^  “  v2t03>  anci  e^'fec't  of  gravity  both  negligible 

aiid  assuming  m  constant  and  6  small.  The  equations  then  become 


*  Fc 

vc  +  icolvc  -  ivl“c  =  5T 


JT  Jr> 

“  - “)v,V 

v  m  *  1  c 


co  +  iku),co  = 
c  1  c 


j Mv, v„ , 
H  1  c 


The  first  and  third  equations  of  this  set  can  easily 
be  solved.  The  solutions  are 


t  -  b>. 


where  the  constant  of  integration  b  is  determined  by 

vo  =  °°th 

and  . 


M1  = 


-r-  dt  +  o>  . 
A  o 


The  second  and  fourth  equations  are  not  easy  to  handle 
because  of  the  variable  coefficients.  It  is  interesting  to 
note,  however,  that  they  can  be  solved  explicitly  by  assuming 
co^  =  av^  and  introducing  the  new  variable  s  defined  by  ds  =  v^dt. 


1 


The  assumption  =  av^  is  perhaps  a  dubious  onej  but  it  is 

not  pointless  to  solve  the  equations  on  the  basis  of  this 
assumption  for  if  a  =  0  is  substituted  in  the  solution  so 
obtained  it  reduces  to  the  solution  for  the  non- spinning  case, 
If  there  =  av^,  and  s  = j' v^dt  are  used,  the  second  and 

fourth  equations  become  ** 


dv  JT  +  Jn  F 

+  (ia  +  ■■  ■  ■  ■ — ii)  v  -  ico  =  — - 

ds  x  m  .  t  c  c  mv. 


*»_  j  j  v  iP o 

.  +  ikaco  -  1JMvc  = 

ds  c  — Bv1 


which  can  be  solved  readily.  The  solutions  are 


v  =  — - —  ■ 

c  Va  2  ’  t 

a0s  f  -cus  _  „ 

-e  J  e  + - £ ^  dt 

Jo  m  ’  ;a  i<  ) 


\rt 

ais  r 

e  )  e 


—  OL-i  S 


ai  iJMFc  ,  *<VJD> 


Bm  +  7  i  “a  +  1 


•  A  ^ 

LalJ  h  : 


a)  =,,~ — 

C 

“2s  f  "“2s  |ij  F0  |i(JL+JD)  V7 

"e  ^ 6  +  (  HH2  -  a  +  ia2j  Hj  « 


in  which 


•Ai-vii 

2 
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where 


^  -  ia(k+l)  ♦  (iliaj 

j32-is.^  +  ika(i4iB). 


The  above  solutions  can  be  used  to  calculate  the  yaw, 
the  angle  of  deviation  and  the  deflection.  For  example,  take 
the  special  case  for  which  =  Oj  F-^  =  cons' t; 

Fc  =  0;  ?^  =  0;  P ^  =  F-pj  =  0.  The  roots  and 


are  then 


al  ~  1 


’M 


B 


arid  hence 


Vc  =  - 


v 

V^B 


sm 


“2  =  -# 

js(t)  -  S (T)]  dl 


03  = 

C 


Fl6 


cx  ,/r 

r— 

[  cos  VJ 

dT. 


Using  the  above  equations,  it  is  easy  to  verify  that 
«2  =  0;  v^  =  0.  The  motion  is  therefore  planar,  and  it  can 

be  assumed  without  loss  df  generality  that  the  rocket  moves 
in  the  c^  c.2  plane. 


The  yaw  is  given  by 

rt 


tan&  = 


vc 

Fle  ( 

V1 

sin 


m 

2F, 


iff 


M  I  Fit  2  F,T  I 

A  '(-i-  +  v  )  -  +  V  ) 

v  m  o'  v  m  o' 


m 


dT. 


This  can  be  written  in  a  simpler  and  more  usable  form  by 
introducing 
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G  ^ 


h  „ .  1/5 

m  ^  1 B 


k2 

*1 


B 

m 


#-l  =/|  (Gt  +  V0)  =  VjTz 


|/§(G  T  +  v0)  =  .  V&*. 


Th^  expression  for  the  yaw  then  becomes 


tan  b  = 


p  k?z 


f 


>in  5(z2-  x2)  dx. 


The  angle  <p  is  given  by  tan  cp  =  jr^/r-J  or  (see  last 
section)  by  tan  cp  =  |  i  J  cocdtj=  |  Vtt/h  G  cocdzj.  Now 


%  =  4^ 
0  k! 


l /^T  f*  cos  §(z=  -  x*)dx. 


An  integration  then  yields 


tan  *  ‘fjfc* 


r,. 

(cos  | 


x2  dx  y  + 


sin  ^  x2  dxf 


'e 

4 


For  a  more  detailed  analysis  of  the  last  results  see  the 
report  by  I.  Bowen,  L.  Davis,  L.  Blitzer:  The  Effect  of  Fin 
Size,  Burning  Time,  and  Projector  Length  on  the  Accuracy  of 
Rockets,  NDRC,  CIT-JPC3. 


Arthur  S.  Peters 
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